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ABSTRACT 

Current mass limits allow the possibility that squarks may be produced in large 
numbers at the next generation of linear e'^e'^ colliders. In this paper we investigate 
the prospects for precision studies of squark masses at such colliders. We assume that 
squarks are lighter than gluinos, and discuss both direct and cascade decay scenarios. 
By exploiting the clean environment and polarizable beams of linear e'^e~ colliders, 
we find that squark mass determinations at the level of a few GeV are possible in a 
large part of the parameter space. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is at present a promising theoretical framework for physics be- 
yond the Standard Model In addition to being free of quadratic divergences and pro- 
viding a natural dark matter candidate, the simplest supersymmetric theory, the Minimal 
Supersymmetric Standard Model (MSSM), has had remarkable success in explaining the 
unification of coupling constants [Q. To solve the gauge hierarchy problem, SUSY must 
be broken at energies of order 1 TeV, and thus the supersymmetric particles of the MSSM 
must be within reach of the next generation of accelerators. This has stimulated a great 
deal of activity in MSSM phenomenology in recent years. It should be noted, however, that 
most of this activity has centered on particle searches. If supersymmetry is discovered, there 
will be a rich spectrum of superparticles, and detailed studies of their masses and couplings 
will be the focus of experimental particle physics into the next century. Precision measure- 
ments of such quantities are crucial to the understanding of the SUSY breaking sector of 
the MSSM and may even inform attempts to understand the mechanism of SUSY breaking 
in supergravity and string theories. It is not, then, premature to investigate the prospects 
for detailed study of superparticle properties at future facilities. 

SUSY events are commonly characterized by unobservable particles in the final state, and 
it is therefore not obvious that the underlying SUSY parameters may be extracted from them 
with enough precision to be of theoretical interest. To date there has been relatively little 
work in this direction, with the exception of recent precision studies of sleptons and neutrali- 
nos/charginos These studies have shown that if such particles are pair-produced 

at a future e+e~ collider, their masses could indeed be determined with enough accuracy 
to be significant tests of grand unified supergravity models. These particles are a natural 
starting point, as they are expected to be the lightest in the supersymmetric spectrum and 
therefore the most likely to be accessible at e+e~ colliders with ^/s = 500 GeV — 1 TeV. 
It is also possible, however, that squarks may be light enough to be produced at such 
machines. This scenario is allowed by current experimental limits and, if true, would 
provide an extremely fertile ground for future experiments. Present limits allow squarks 
with masses above 100 GeV. With the expected features of the next e"'"e~ collider, luminosi- 
ties of 10fb~^/year and beam energies of 250 GeV [3,p|,p!0|, these would be pair-produced 
at the rate of thousands per year. More massive squarks would have to wait for accelerator 
upgrades, but their study would be qualitatively similar. Although squark searches have 
been considered in great detail in the context of hadron colliders [lTTl , p!2| , p!3| , |l^ , [T5| , p!6| , pre- 
cision studies at these machines would be very difficult. Not only must one control large 
hadronic backgrounds, but one must also work in a situation where only a few features of 
a complex event are reconstructed. In contrast, the clean environment of e~^e~ accelerators 
makes them promising for such studies. An added advantage is the availability of polarized 
beams, which, as we will see, is a very useful tool. The purpose of this article is to survey 
the prospects for precise determination of squark masses in the e^e~ environment. 

The discussion is organized as follows: In Sec. |T| we review relevant information about 
squarks and the MSSM. We will note precisely which simplifying assumptions of the MSSM 
we use in our analysis. In Sec. |Tl| we describe our event simulation and organize our 
exploration of the SUSY parameter space. In Sec. ^ we treat the simplest region, in which 
squarks decay directly to the lightest supersymmetric particle (LSP) . The more complicated 
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regions with cascade decays are investigated in Sees. and 0, where the hght neutrahnos 
are gaugino-hke and higgsino-hke, respectively. In Sec. [VII| we extend the analysis to the case 
of higher mass squarks and higher energy colliders. We conclude with some final remarks in 
Sec. WTB. 



II. SQUARKS IN THE MSSM 

The MSSM is the simplest extension of the Standard Model that incorporates super- 
symmetry. Here we concentrate on the salient features for squarks and their decays. Full 
discussions may be found in a number of articles . 



We will use hats to denote superfields and tildes to denote superpartners of Standard 
Model particles. The MSSM includes the usual matter superfields and two Higgs doublet 
superfields 

where Hi and H2 give masses to the isospin — | and +| fields, respectively. These two 
superfields are coupled in the superpotential through the term —fieijHlH^, where /i is the 
supersymmetric Higgs mass parameter. Both /i and tan/3 = {H2) / {Hf) will be important 
for our analysis of squark decays. Soft supersymmetry breaking terms for scalars and 
gauginos are included in the MSSM with 



3 



l^soft = T.M^\<P^? + \{ WiBB + ^ M^Wm^ + M^-g'^t] + h.c. ) , (2) 

^ ' 0=1 k=l 



where i runs over all scalars. 

For every fiavor of quark in the Standard Model, the MSSM contains both a left- and a 
right-handed squark whose masses are given by 

= M? + + m| cos 2(3 {\- fsin^^iy) 
m?^ = M? +m2 + |m|cos2/?sin2^iy 

= M^ + ml-ml cos 2/3 (| - Isin^Ow) 
= M|j + — |m| cos 2(3 sin^Ow, 

where we have suppressed generational indices and have neglected left-right mixing terms 
which are relevant only for the top and bottom squarks. The terms Mq,Mq, and are 
some of the soft supersymmetry breaking terms Mj of Eq. |^, and the remaining terms are 
model independent and known once tan j3 has been determined. Accurate determination of 
squark masses thus allows determination of the SUSY breaking parameters of the squark 
sector, which are of great theoretical interest. 

In the most general form of the MSSM, the different SUSY breaking parameters are un- 
related. To satisfy fiavor changing neutral current constraints, however, corresponding terms 
in the first and second generations must be nearly degenerate (Recent work has noted 



that adequate squark degeneracy may be enforced by gauged horizontal symmetries pO| or 
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may in fact be unnecessary to satisfy the FCNC constraints [^.) In minimal low energy 
supergravity models, the SUSY breaking parameters are usually related by assuming that 
they evolve from a universal scalar mass at a high scale. Renormalization group evolution 
then predicts a splitting between the left- and right-handed terms due to the contribution of 
SU (2) gaugino loops, which can be of order 5% or higher depending on the gaugino mass, 
as well as a splitting of the third generation squarks from the first two. However, this as- 
sumption of a universal scalar mass is not on firm theoretical footing, and one of the goals of 
precision mass determinations is to differentiate between the universal scalar mass scenario 
and one where the different masses are not so simply related. We therefore do not assume 
a universal scalar mass. 

The experimental signature of squark production will depend crucially on the details of 
the squark decay channels. As we will see below, the particles of the Higgs and gaugino sec- 
tors appear as intermediate decay states, and it is therefore necessary to discuss these sectors 
of the theory. The Higgs sector consists of two CP-even scalars and the charged scalar 
H^, the CP-odd scalar and the Higgsinos iJf, H^i and iJg- The Higgsinos mix 
with gauginos to form mass eigenstates. The charged mass terms are (?/'~)^M£±'?/'"^ -|- h.c, 
where {il)^Y' = {—iW^,H^) and 

f M2 V2Mw sin f3\ 

The chargino mass eigenstates are xt = ^iji^t Xi~ = Uj-,?/'", where the unitary ma- 
trices U and V are chosen to diagonalize M-j^±. Neutral mass terms may be written as 
i(V'°)^M^oV° + h.c, where = i'^B, -iW^ H^, H^) and 



I Ml -Mzcos^sin6'vy MzsmpsmOw \ 

_ M2 cos /5 cos 6'vi/ — sin/5cos6'vy 

^ —Mz cos [3 sin 6^ Mz cos /5 cos 9]^ — 

V sin /? sin 6'iy — sin/3cos6'vi/ — / 



(5) 



The neutralino mass eigenstates are = Nijipj, where N diagonalizes M^o. The four 
neutralinos in increasing order of mass are Xi^ X25 X35 X% ^iiid the two charginos, 
similarly ordered, are xt and xt- 

The MSSM as reviewed above contains many arbitrary constants, and usually some sim- 
plifications are made. In our version of the MSSM we will employ the following assumptions: 

(a) We assume R-parity conservation, a common assumption that imposes baryon and lepton 
number conservation. R-parity conservation implies that the LSP is stable and must be 
among the decay products of any supersymmetric particle. 

(b) We take the LSP to be the lightest neutralino, Xi- The LSP must be neutral, and the 
other candidate, the sneutrino z/, is heavily disfavored if one assumes that the LSP makes 
up galactic dark matter |2^. With this assumption, squarks can either decay directly into 



an LSP and a quark, 

Q Qxl, (6) 
or indirectly to an LSP through a chain of neutralinos and charginos, e.g., 
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q q'xt Q'XiiQQ, ^l) or qxl qXiill, ^^)- (7) 

The latter decays are called cascade decays. Even more complicated cascades involving x% 
X4, and X2 are possible. The LSP interacts very weakly and disappears from the detector 
like a neutrino, so direct LSP decays leave a distinctive signature of acoplanar jets + missing 
Pt(/t)- In the case of cascade decays, there may be additional jets and leptons. 

(c) We assume rriq < IM3I so that squarks do not decay through gluinos. Without this 
assumption, the decay q ^ qg would be possible and would in fact dominate. Since the 
gluino is strongly interacting, it would be accompanied by additional hadronic radiation and 
a weakened spectrum. In such different analysis from ours will be needed. 

(d) In this discussion, we assume that the neutralino/chargino sector has already been 
explored. We will therefore keep this sector simple by assuming the unification of gauge 
constants and gaugino masses at some higher scale. This implies that even at lower scales 
we have the one- loop condition |^| 



M2 ^ 3 Ml ^ M3 

9l 5 g'^ gl ^ ^ 

Eq. D, coupled with assumption (c), will make our analysis invalid in some areas of the 
(/i, M2) plane, which will be noted below. 

Because direct decays of the squarks to the LSP are characterized by two acoplanar jets 
with large /t, they should be easy to resolve experimentally. For example, cutting events 
with < 35 GeV and 9acop > 150° will eliminate the bulk of Standard Model backgrounds 
23] while eliminating only 20% - 40% of the signal events. The primary backgrounds after 



these cuts are W^W~ production, e^pW"^ via jW fusion, and uuZ via WW fusion. W pair 
production is a background to direct squark decays when one W decays leptonically and 
the charged lepton is either mistakenly included in a jet or goes undetected. In the former 
case one can eliminate events whose missing momentum and energy are consistent with an 
undetected neutrino, and in the latter case one can eliminate events whose visible mass is 
consistent with the W mass. All in all, we expect W pair production to be a less troublesome 
background for squark studies than it is for chargino and slepton studies. A more serious 
background is e^uW"^. Before cuts, this cross-section is roughly an order of magnitude above 
the signal, and because the electron tends to be lost down the beam pipe and the neutrino 
tends to be produced with large pr, the cuts are not particularly effective. However, this 
background can be removed by a cut on the two jet invariant mass. By eliminating events 
with invariant mass less than 100 GeV the background is effectively removed while much of 
the signal is retained. Such a cut also reduces the uuZ background to negligible levels. One 
might also be able to take advantage of the fact that these backgrounds tend to produce jets 
in the same hemisphere, while squark pair jets are preferentially in opposite hemispheres. 

Cascade decays result in a wide variety of signals, and are therefore generally more 
difficult to isolate. However, in the analysis of this paper, we will concentrate primarily 
on events with two jets and isolated leptons. Backgrounds to such events consist largely 
of events where an on-shell W or Z decays hadronically, and these will therefore also be 
removed by a cut on the two jet invariant mass. One important exception to this is ti 
events, which should be well-understood. This is a background when both top quarks decay 
to blu. We can reduce this background by anti-tagging b quarks, and in our discussion of 
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cascade decays, we will discuss the effectiveness of such a cut. This cut excludes b and t 
squarks from our analysis: however, for reasons outlined in the following section, the third 
generation of squarks will most likely require a separate analysis anyway. 

In addition to the Standard Model backgrounds discussed above, there may also be SUSY 
backgrounds to consider, such as chargino and neutralino production. Chargino pairs are a 
background to 2-jet cascade events when one chargino decays through a hadronic W decay 
and the other through a leptonic W decay. Neutralino production can result in backgrounds 
to both direct and cascade squark decays. However, in both chargino and neutralino 2-jet 
events, the two jets are produced by an on- or off-shell W and Z decay, so the two jet 
invariant mass should be less than or of order the W and Z masses. Thus, the invariant 
mass cuts that reduce the Standard Model backgrounds should also effectively reduce the 
SUSY backgrounds. In any case, it is not unreasonable to hope that once the Standard 
Model and MSSM backgrounds are well understood, efficient cuts can be devised to isolate 
the squark signal. 

In the remainder of this paper we will apply only the and acoplanarity cuts given 
above, namely ^/r < 35 GeV for all events and Oacop > 150° for events with only two jets. A 
cut requiring the 2-jet invariant mass to be greater than 100 GeV has only a small effect on 
the distributions we will plot, though it decreases the total signal by roughly a factor of 3. 

III. EVENT SIMULATION AND ORGANIZATION OF PARAMETER SPACE 

For the purposes of this exploratory study, we use a simple parton level Monte Carlo 
to simulate squark production and decay. We then simulate hadronization and detector 
effects by smearing quark jets with a detector resolution of a^'^/E = 50%/\/E' {E in GeV). 
The cuts on j/t and 9acop discussed above are implemented to separate Standard Model 
backgrounds. Beamstrahlung and initial state radiation are not included, and we assume 
100% electron beam polarization. 

Systematic errors in this study are of two kinds. There are, of course, systematic errors 
arising from hadronization and detector effects. As our modeling of these effects is rather 
crude, a detailed study of these errors will not be attempted here. In addition, there are 
errors arising from the uncertainty in the determination of the neutralino and chargino 
masses which enter our analysis. The effects of these errors on our ability to determine 
squark masses will be included in our discussion of the case of direct decays. 

We must fix the squark masses for the Monte-Carlo simulations. To simplify the analysis, 
we assume that all left-handed squarks from the first two generations are degenerate, as are 
the right-handed squarks, and concentrate on the left-right splitting. This assumption may 
be relaxed without major changes to the analysis, but will result in greater complications 
in fitting the data. The third generation of squarks will require a separate analysis, as the 
heavy top mass may not allow direct decays of the stop to the LSP, and in addition there 
will be an appreciable left-right mixing. In most situations, the third generation of squarks 
can be separated with b anti-tagging, so we will consider only the first two generations from 
here on. 

The present lower mass limit on squarks is ~ 100 GeV from the CDF experiment [^] when 
the important effects of cascade decays have been included [|1^]. These limits assume M3 < 
400 GeV and disappear for higher M3. Squarks will be pair-produced by e~^e~ —>■ (7, Z) — »• 
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qq. At e~^e~ colliders with y/s = 500 GeV, the production rate becomes substantial not far 
below the kinematic limit of rriq = 250 GeV. For our studies, we will choose m^^ = 220 GeV 
and rrigj^ =210 GeV. These values are significantly above the present mass bounds, but are 
also low enough to give ~ 2500 events per year, assuming two degenerate squark generations 
and unpolarized beams at a luminosity of lOfb^^/year. The results for uig = 400 GeV and 
y/s = 1 TeV are qualitatively similar; we will address this case in Sec. [VII . 

In Fig. |I|, we plot the cross sections for pair production of g^g^ and qnqji from e'^e~ 
annihilation, assuming polarized e~ beams. Note the slow rise of the cross section near 
threshold, characteristic of scalar particle production, which precludes the use of cross section 
measurements alone for precision mass determinations. The polarization dependence of the 
cross section is an important feature. We see that beams produce qiqi pairs 91% of the 
time, and beams produce qRqpi pairs 91% of the time. Beam polarization is therefore an 
extremely effective way of separating the left- and right-handed squarks. 

Since squarks are expected to be among the heaviest superparticles, one cannot ignore 
the possibility of cascade decays. We should therefore consider a representative variety of 
values of the neutralino/chargino sector parameters rather than just the limiting cases. To 
do this, we must organize our survey of the parameter space. Given the constraint (H), all 
chargino and neutralino masses are given as functions of tan/5, /i, and one of the gaugino 
masses, which we will take to be M2. Thus, for a given value of tan/5, squark decays are 
determined by the parameters (/i, M2), and our task reduces to exploring squark decays in 
this parameter plane. 

The discussion will be limited to the part of the plane with < M2 < 1 TeV and 
|/i| < 1 TeV. The constraint M2 > may be imposed without loss of generality. The upper 
bound on M2 results from the fact that M2 is a SUSY breaking parameter, and for SUSY to 
naturally explain the electroweak scale, M2 cannot be too large. Although fi is not a SUSY 
breaking parameter, a large value for fi would re-introduce the fine-tuning problem. For 
most of this study, we will fix tan (3 = 2. The analysis for higher tan (3 is not much different 



and will be deferred to Sec. VIII 



With the parameters chosen above, we may now begin our study in earnest by dividing 
the (/X, M2) plane into regions with similar squark decay channels. The decay patterns 
will, of course, also be influenced by the Higgs and slepton masses. Since these masses are 
unknown, however, the boundaries separating regions with different squark decays through 
Higgs and sleptons are not fixed. For this reason, we first divide the parameter space into 
regions ignoring the Higgs and slepton decays; we will then consider the effects of varying 
the Higgs and slepton masses within each region. The regions are shown in Fig. ^ for 
rriq = 220 GeV. We can ignore the hatched regions in the upper left and right corners, 
because there the squarks are lighter than x^, which violates our assumptions. The cross- 
hatched region along the M2 = and /i = axes is experimentally ruled out by lower 



bounds on superparticle masses p5| , |26| . In region 4, the squarks can decay to four or more 
of the neutralinos/charginos. This results in very complicated cascade patterns. Though 
our methods may, in principle, be extended to this region, we will not consider it further 
in this paper. In what remains of the plane, the situation is greatly simplified, because the 
three neutralinos/charginos to which the squark can decay are always Xi, X2) and xf- The 
only two-body decays kinematically allowed in region 1 are those directly to the LSP. In 
regions 2 and 3, squarks may also decay to either X2 xt or both. In region 2, x1, X2, 
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and Xi cill dominated by their gaugino components. This region is further subdivided 
into region 2a, where X2 Xi may decay to the LSP through on-shell W and Z bosons, 
and region 2b, where decays through on-shell W and Z decays are not possible. In region 
3) Xi; X2; ^^cl xf are all Higgsino-like, and again only decays through off-shell Ws and Zs 
are possible. Finally, each region has a mirror region in the /i > part of the plane. These 
mirror regions usually have the same decay patterns as their /i < counterparts, so where 
this is true, we will analyze only the fi < regions. The only exception is region 3, and we 
will note the different behavior for > at the end of Sec. 0. In Fig. |^ the representative 
points in each region that we will consider in detail in the following sections are marked. 

In what follows, it will often be helpful to keep in mind the following approximate 
relationships [|^. At Mz, Mi ^ ^Ma and M3 ^ ^Ms. For M2 > M^, 

m^o ^ min{|/i|, iMa}, m^o ^ min{|/i|, M2}, 

m^o ^ max{|/i|, IM2}, m^o ^ max{|/i|, M2}, (9) 

m^± ~ min{|/i|, M2}, m^± ~ max{|/i|, M2}. 

Note that in this approximation, X2 ^^^1 xf virtually degenerate throughout the plane. 
We will take advantage of this fact in our analysis of cascade decays in Sec. 0. 

It should also be noted that our assumption that the squarks are lighter than gluinos, 
coupled with the unification assumption of Eq. ^ implies that our analysis is only valid for 
M2 > 67GeV (above the dotted line in Fig. |). Also, if one makes the further assumption 
that one can compute the squark mass by applying the renormalization group equations 
with the desert hypothesis, one obtains 

m? ^ml + 7m^/2, (10) 

where mo and mi/2 are the squark mass and the gaugino mass at the unification scale Mu, 
respectively. Since ml > and mi/2 ~ O.8M2, this leads to the constraint M2 ^ O.Sm^, 
which in the present case implies that M2 ^110 GeV. 



IV. DIRECT DECAYS 



In region 1 only direct decays of the squarks to the LSP are allowed. For the sake 
of concreteness, we will perform our Monte-Carlo simulations at the representative point 
(/X, M2) = (—500 GeV, 300 GeV). The quark jet from the decay of the scalar particle q 
will have a flat energy distribution with endpoints 



qxi 




X? 

2 

mi 



(11) 



where is the beam energy, and we have neglected the quark mass. Thus, in theory, 
squark masses can simply be deduced from the distribution's endpoints. Unfortunately, the 
simple fiat shape will be changed by cuts, and finite detector resolution and hadronization 
will smear the endpoints. Of the 1764 (975) squark pair events produced by the (e^) 
beam, 1294 (683) survive the cuts. The energy distribution of the individual jets from the 
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surviving events is given in Fig. where detector resolution effects have been included. We 
see that the location of the endpoints is rather ambiguous. One way to extract the squark 
mass is to apply a binned likelihood fit to the jet energy distribution, with the logarithm of 
the likelihood given by 

ln£(mg^,mg^) = ^ v4j(mg^, m^^) In^j - Bi, (12) 

i 

where the sum is over all bins, Ai{mqj^, niq^ is the expected number of events in bin i given 
hypothetical squark masses rriq^ and m^^, and is the measured number of events in bin 
i. The actual values of rriq^ and determined by maximizing ln£, and the statistical 

error of the determination is given by the width of the ln£ peak. For simplicity we leave 
TTLqj^ fixed at its actual value and calculate ln£(mq^). In the ln£ calculations, we have 
approximated the theoretically expected number of events Ai by Monte-Carlo simulations 
with a very large number of events (typically 50,000). We find that the m^^ determination 
has a statistical error of 0.9 GeV at 95% CL. These calculations have been performed 
assuming that the LSP mass is known. The expected statistical errors on the masses of 
xf and Xi from chargino studies are 3.2 and 2.0 GeV, respectively, while slepton studies 
should give the LSP mass to 1.0 GeV A shift of 1.0 GeV in the LSP mass causes a shift 
of about 0.7 GeV in the central value of the likelihood fit. 

In the single jet energy spectrum used above, correlations between the energies and 
directions of the two jets of a given event are ignored, since each point in the distribution 
represents only one jet. It is thus worth thinking about whether there is a more efficient 
way to use the information contained in the event sample. One possibility is to retain the 
jet correlation information by using a two-dimensional binning, but this necessitates a large 
jump in computing time to obtain the distribution accurately. 

Instead, the method we will use extensively is the following: For each event, we calculate 
the quantity m™™, the minimum squark mass kinematically possible, given the two observed 
quark jet momenta. For each event, we then get one value of m™™, and we apply a binned 
likelihood method to the m™™ distribution. The quantity m?™ is easily determined (see 
Fig. H). We label the particle momenta as q{pi) Xi(P3)Q'(P4) and qip'i) — > Xi(P5)9(P6)- 
The total visible momentum is then py = p4 +P6- Because we know and have measured 
E\ and £^2, we can determine the LSP energies £"3 and E^. However, we also know the LSP 
mass, so we can find the magnitudes \'pz\ and jpsl, and since ps -|- = — py, the vectors 
P3 and ps are constrained to lie on a circle C . We then can calculate the angles 7 and 5 
shown in Fig. ^. The minimum squark mass corresponds to the maximum possible |p2|, and 
is given by 

{mf^f = El - |p3|2 - |p4|2 + 2|p3||p4|(cos7cos5 - sin7sin5). (13) 

The distributions of m?™ for left- and right-polarized electron beams are sharply peaked at 
the underlying masses rUq^ and m^^, respectively (see Fig. ||). As before, for simplicity we 
leave fixed at its actual value and calculate In C{mq^) . We find that the statistical error 
of the rriq^ determination is reduced to 0.4 GeV at 95% CL. Thus, we see that the m™™ 
method improves our statistical error significantly. If the cut m2_jct > 100 GeV is necessary 
to reduce background, our event sample is reduced by a factor of 3, but the statistical error 
increases only slightly to 0.5 GeV at 95% CL. We note, however, that a shift of 1 GeV in 
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the LSP mass now shifts the central value of the likelihood fit by approximately 1.5 GeV. 
The m™"^ calculation is more strongly dependent on the LSP mass, and therefore places a 
higher premium on its accurate determination. 

The greater power of the m™''^ distribution is manifested in its sharp peak, because even 
slight variations in m^^ and move these peaks enough to create large differences between 
Ai and Bi in some bins. This implies that \nC falls rapidly from its maximum. There are 
a number of reasons for the sharp peak of the m™™ distribution. Here we just note that, 
roughly speaking, momentum vectors lying on large circles C may give mass minima both 
close and far from the actual squark mass, depending on where the momentum vectors lie 
on C. However, small circles give only accurate solutions, and thus the calculated minimum 
masses preferentially lie close to the actual underlying squark mass. 

One potentially powerful feature of squark mass studies in this region is that, since 
both left- and right-handed squarks have identical decay channels, a direct comparison can 
be made to determine left-right mass splittings. The left- and right-handed squarks can be 
isolated using polarized beams, and systematic errors, which should effect both polarizations 
equally, should largely cancel in the ratio of their masses. One can therefore determine left- 
right mass splittings to greater accuracy than one can determine the actual values of the 
masses. This is in contrast to the case of slepton studies, where use of the left-polarized 
beam is hampered by a large W~^W~ background explained above in Sec. 0, 

we do not expect W^W~ backgrounds to be a signifigant problem for squarks. 

While it may come as no surprise that precise mass determinations can be made in the 
simple case of direct decays, the large squark masses expected make it likely that more 
complicated decays will be present. For a generic point in parameter space, with a large 
number of possible decay chains present simultaneously, it is important to determine whether 
it is still be possible to extract accurate squark masses from the more complicated signal. 
We now turn our attention to this question in representive regions of parameter space. 



V. GAUGINO-LIKE CASCADES 

In region 2 squarks have new decay channels through on-shell X2 ^i-nd Xi- this region, 
however, B, y^K, W^, and xf ~ W^, so since right-handed squarks do not couple to 
SU{2) gauginos, they still decay predominantly directly to the LSP. This may be seen from 
the branching ratios of ul and £t/j, given in the contour plots of Fig. ^. Similar plots for 
down-type squarks differ little. The analysis for right-handed squarks is therefore simple. 
First we use the polarized beam to isolate right-handed squarks. Contamination from 
left-handed squarks will be of order 10%, and most of these will go through complicated 
decay channels and can be easily separated by considering only two-jet events. We then 
apply the analysis of region 1 with little degradation of statistics. As shown in Fig. ^, left- 
handed squarks, unlike right-handed squarks, do decay predominantly through cascades in 
some parts of region 2. When this is the case, a separate analysis for left-handed squarks is 
necessary. The rest of this section will be concerned with left-handed squarks only. 

As suggested at the end of Sec. |T|, we can use the near-degeneracy of X2 xf to 
simplify matters. Throughout region 2, these two particles are typically degenerate to a few 
GeV, and since 
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phase space suppression allows us to safely ignore decays of X2 xf ^-iid vice versa. We are 
then left with two-step decays with squarks decaying to X2 ^^nd xf, which then decay to the 
LSP through three-body modes mediated hj W or Z bosons, sleptons, sneutrinos, squarks, 
or Higgs bosons. The resulting quark jet energy distribution is much more complicated than 
it was in the direct case, with quarks produced both at the primary vertices, i.e., the initial 
squark decay vertices, and in the later cascade decays. 

However, by choosing appropriate cuts, we can reduce the problem to the case of direct 
decays. Recall that in the m?™ analysis of region 1, we needed to know only the energies 
and momenta of the quarks leaving the primary vertices and the mass of the neutralino or 
chargino leaving those vertices, which we now denote m^^^'^'^^ . In the case of region 1, such 
identifications are obvious, since in all events both of the quark jets are produced at the 
primary (and only) vertices, and m?"™^'^'^ = m^o is always the correct choice. Cascade decays 
complicate this analysis, because the initial squark decay may involve neutralinos/charginos 
other than the LSP, and there may be many quarks produced, making it difficult to determine 
which two came from the primary vertices. Our stategy will be to find cuts that allow us to 
accurately identify the primary quarks and assign m^^'^^^^ . We can then calculate m™™ for 
each event and proceed as in the previous section. Two simple strategies are (1) separating 
direct and cascade decays kinematically and choosing events that contain direct decays on 
both sides, and (2) considering double cascade events in which the X2 xf has only leptonic 
decay products. We will show that these two strategies are always sufficiently effective. We 
now consider regions 2a and 2b in turn. 

Decays of X2 ^^'^ xf through on-shell Ws and Zs are allowed in region 2a. We will take 
the representative point to be (/x, M2) = (—500 GeV, 200 GeV). At this point, the masses of 
x1, X2, and xf are 103.2, 206.1, and 206.0 GeV, respectively. The jet energy distribution 
is shown in Fig. |^, where every quark jet produced in the Monte-Carlo squark decays is 
binned separately. The distribution consists of three parts — a large peak of soft quarks 
emanating from squarks decaying to X2 xf (long dashes), a flat distribution from direct 
LSP decays (solid), and a wide hump from hadronic W and Z decays (short dashes). The 
total distribution is given by the dotted histogram. The soft jets in the peak may not be 
discernible experimentally, but our analysis will use only jets from direct decays, and these 
jets have energies that are always greater that 40 GeV. 

As noted previously, since region 2 allows cascade decays, top quark production is poten- 
tially a troublesome background. In Fig. ^, we again plot the total jet energy distribution 
resulting from squark decays, but this time, for comparison, also the bottom quark spectrum 
resulting from top decays. In the bottom quark spectrum, we have removed all b quarks 
produced in top events in which at least one b quark has been successfully tagged. We 
assume mtop = 150 GeV and a 6-tagging efficiency of 80% 0. The top quark background 
is substantially reduced below the signal, and as top quark decays will be well- understood 
by the time of these squark studies, the top quark background should not be a significant 
obstacle. 

The branching ratios at (/x, M2) = (-500 GeV, 200 GeV) are 
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BR{uL ux?) = 58%, BR{dL ^x?) = 36%, 



BR{ul ux^) = 28%, BR(4 dxs) = 43%, (15) 
BR{ul dxt) = 14%, BR(4 ^ uxr) = 21%, 

so we see that there are still many direct LSP decays. We will therefore try to isolate the 
double direct LSP decay events and then apply the analysis of Sec. fV[ Of the 1764 events 
produced with a left-polarized beam in 1/2 year, 1437 pass the and 9acop cuts. We then 
consider only events with 2 jets and no isolated leptons, leaving 673 events, and for each of 
these we calculate m?™ using m|™^'^'^ = fn^o. Of the remaining events, 244 are actually 
cascade events with neutrino decay products, but these may be removed by considering only 
events in which both jets have energy above 30 GeV. We find that the likelihood fit to the 
remaining m™™ distribution gives squark masses to 1.3 GeV at 95% CL. 

The distribution of Fig. |^ was calculated for the case where the sleptons and Higgs scalars 
are massive enough that their diagrams are off-shell and suppressed. Smaller slepton and 
Higgs masses will, of course, change the part of the jet distribution resulting from cascades, 
but will have no effect on the two-jet event sample and will therefore not change our results. 

We now turn to region 2b, whose representative point we take to be (/i, M2) = 
(-500 GeV, 100 GeV). Here the masses of x?, X2, and xf are 52.8, 108.1, and 107.8 GeV, 
respectively, and we see that the W and Z diagrams are indeed off-shell. The branching 
ratios are 

BRiuL^uxD = 7%, BR{dL ^ dxl) = 1%, 

BR{ul ^ MX2) = 61%, BR{dL dxl) = 67%, (16) 

BR{ul dxt) = 32%, BR{dL ^ uxi) = 32%. 

Because there are not many direct LSP decays, the analysis must rely on cascade events. 
This complicates the analysis. However, we now show that even for the most difficult sets 
of parameters, a significant number of X2 and xf decays will be purely leptonic. In these 
events we may unambiguously identify the two primary vertex quark jets, and we can then 
again apply the region 1 analysis. 

To do this we must analyze the relative importance of the various cascade diagrams (see 
Fig. ^). The Higgs diagrams may be safely ignored. Typical values for the mass of the 
lightest Higgs scalar, h'^, are in the range 70-110 GeV. For this entire range, the process 
proceeds off-shell and is also suppressed by the bottom quark Yukawa coupling. The other 
Higgs scalars are significantly more massive and can also be safely ignored. In addition, we 
have squark and slepton cascade diagrams. Sleptons are generally expected to be lighter 
than squarks, and we will see below that the lower the slepton masses, the simpler our 
analysis. For now, we will pessimistically take the sleptons to be degenerate and of mass 
200 GeV. 

Since we have squark and slepton masses significantly higher than M^, one might expect 
the W and Z diagrams to dominate. However, in the gaugino-like region of parameter space, 
|Nii| ^ |Vii| ^ 1, so the LSP is primarily composed of B, which does not couple to the W or 
Z at all. Thus, these diagrams are suppressed by a factor S = 0(|Vi2|, IN12I, IN13I, |Ni4|). 
We can obtain a rough estimate of 5* by taking an appropriate limit of the explicit expression 



for V12 |2^. Taking Mw small compared to /i and M2, we have 
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for tan (3 = 2, where we have expanded in terms of M2//U, a small parameter in the gaugino- 
like region. Thus, V12 is about 0.02 at our representative point. The Ni^ are roughly of the 
same order, and we find that 0.01 ^ ^ 0.1 throughout most of the region. The competing 
suppressions of the different diagrams make it impossible to simply determine which diagram 
dominates. In fact, our calculations show the gauge boson, squark, and slepton diagrams to 
be roughly of the same order. (Note that if one takes the gaugino-like relations Xi ~ —iB, 
X2 ~ —iW^, and xf ~ —iW^ too seriously, one is led to incorrectly conclude that the W 
and Z diagrams may be set to zero. This is never valid in the part of the plane we are 
considering.) 

Evaluating the various diagrams with mj = iriy = 200 GeV, we find 

BR(x^ ggX?) = 20%, BR{xf -> g'gx?) = 55%, 

BR(x^ ^ 11x1) = 31%, BRixt ^ li^x'i) = 45%, (18) 

BR(x^ ^ z/z/X?) = 49%. 

We see that there are many leptonic decays, and we can therefore consider only the events 
with two quark jets + leptons without a great loss in statistics. (Note that at this point in 
the parameter plane the only jets below 40 GeV are those produced by Ws and Zs. Thus 
by selecting only two-jet events, the jets that enter our analysis are again all hard enough 
to be experimentally detected.) Of the 1508 polarized events that pass the ]/t and 9acop 
cuts, 834 have exactly 2 quark jets (+ leptons). Because the polarized beam has eliminated 
most right-handed squarks from our sample, and left-handed squarks rarely decay directly 
to the LSP, we expect that most of these events are double cascade events, and in fact 
70% are. We then calculate m™" with the assignment m?'^™'^'^-^ = 'm^±. Note that since X2 
and xf are virtually degenerate, we need not distinguish them for the kinematic analysis 
of Sec. 1^. About half of the events with one or two direct LSP decays are kinematically 
incompatible with the assignment m^"'^'^^^ = m-±, and we eliminate these. Almost none 

A. X.\ 

of the double cascade events are removed by this cut. We are then left with 689 events, 
of which 84% are double cascade events. The double cascade events form a sharp peak in 
the m™™ distribution, while the few events involving direct decays are much more broadly 
distributed. After performing a likelihood fit to the m?™ distribution, we find that the 
squark mass can be determined to 2.4 GeV at 95% CL. 

If we lower the slepton mass, the slepton diagram contribution grows. For mj =110 GeV, 
even though the sleptons are still off-shell, the gauge boson suppression factor S defined 
earlier allows the slepton diagrams to dominate, and the resulting branching ratios are 

BR(x^ ^ qqx'i) = 0.3%, BR(xf ^ q'qx'i) = 6%, 

BR(x^ ^ /fx?) = 39.8%, BR(xf ^ lux^i) = 94%, (19) 

BR(x^ ^ uux'i) = 59.9%. 

Clearly most X2 and xf decay to leptons, and we can use almost all of the cascade events 
in our analysis. 

VI. HIGGSINO-LIKE CASCADES 

In region 3, xt xf are again close in mass, but now, as suggested in Eq. |, they 
are also close to the LSP in mass. We will take as our representative point (/x, M2) = 
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(-100GeV,700GeV), where the masses of x?, X2, and xf are 98.3, 110.9, and 106.1 GeV, 
respectively. The branching ratios of left- and right-handed squarks to the LSP are now 
both ~ 10% - 20%. The analysis for left- and right-handed squarks is similar, and we will 
consider only the left-handed below. 



The quark jet energy distribution is shown in Fig. The large low energy hump is 
composed of the soft quarks produced in the decays between the neutralinos and chargino. 
The primary vertex quarks have the flat energy distribution we expect, but this is really a 
superposition of decays to three different particles with slightly different masses. We see that 
the primary and secondary vertex quark jets are well-separated in energy, and we can isolate 
the primary quark jets with a simple jet energy cut at 30 GeV. In fact, we again have not 
included the decays of X2 Xi the Monte-Carlo simulation. This omission is physically 
unwarranted, as such decays no longer suffer the pronounced phase space suppression relative 
to the decays of X2 and xf to the LSP. However, these decays will only increase the number 
of soft jets in the low energy hump, which will be eliminated with the energy cut anyway. 
We have also assumed m^- = 200 GeV and m^o = 110 GeV. Other values will change the 
shape of the low energy hump, but again this is irrelevant after the energy cut. 

We must now determine the squark mass from the jet energies. Unlike in the previous 
sections, we cannot use the m™™ distribution, because for a given quark jet, we cannot tell 
if its primary vertex partner was X2 0^ xf^ and the 5 GeV mass difference between these 
two is now significant compared to the the accuracy with which we hope to measure the 
squark masses. We will therefore simply use the jet energy distribution for our likelihood 
fit. Our strategy will then rely on the assumption that we know the neutralino and chargino 
parameters. The x° and x^ masses will be quite accurately determined and will fix the 
relative positions of the endpoints of the energy distributions. Slepton and chargino studies 
should be able to determine the parameters /i and M2, which will determine the branching 
ratios of the squarks. In the Higgsino-like region, it may be difficult to determine M2 
precisely, but since the branching ratios are relatively insensitive to M2, we nevertheless 
also expect to determine the branching ratios accurately. Given these assumptions, we take 
the neutralino/chargino masses and branching ratios as inputs in our analysis. The only 
unknown is then the squark mass, and a In £ fit gives us the squark mass to within 1.2 GeV 
at 95% CL. 

Region 3 is the only region where the corresponding yU > region has different decay 
properties. In the /i > region, all types of squarks decay to the LSP with branching ratios 
^ 80%. We can therefore use the direct decays in this region, and the analysis is actually 
simpler. 



VII. HIGHER SQUARK MASSES 

So far we have been studying squarks with mass ~ 220 GeV. Of course, squarks may 
be significantly more massive than this, and it is therefore important to consider the ap- 
plicability of the preceding analysis to the case of higher mass squarks and higher energy 
colliders. With slight modifications, we will see that it is straightforward to adapt the anal- 
ysis to higher energies. We note, however, that as the squark mass rises, our assumption of 
■rriq < M3 becomes more and more disfavored if one subscribes to the theoretical prejudice 
of a universal scalar mass. 
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To proceed, we will consider the case of a 1 TeV collider with luminosity 30fb~^/year. 
Squarks with mass near the kinematic limit of 500 GeV can be studied; we will take rriq^ = 
400 GeV and m^^ = 390 GeV. The (/i, M2) plane may be divided into regions as before 
(see Fig. |TT]). All of the region boundaries move to higher and M2, except the boundary 
between regions 2a and 2b. The boundary at which gluinos become less massive than squarks 
moves up to M2 ~ 121 GeV. In addition, renormalization group equatons coupled with the 
desert hypothesis now imply M2 ^ 200 GeV. Of course, region 4, the region of complicated 
cascade decays, becomes large. A full treatment of this region would demand a significantly 
more complicated analysis. 

The squark decays of regions 1, 2b, and 3 are qualitatively similar to those in the lower 
energy case we considered earlier, and therefore similar analyses are applicable. An impor- 
tant difference for our analysis, however, is that region 2a is no longer a thin strip, and 
consequently the branching ratio to the LSP is small throughout most of region 2a. In 
the previous discussion, the branching ratio to the LSP, even for left-handed squarks, was 
substantial in region 2a, and we based our analysis on the abundance of direct LSP decay. 
We see that this convenient feature does not persist to cases of higher squark masses and 
higher energy colliders, and for the case of 400 GeV squarks we must use cascade decays for 
left-handed squarks in region 2a. 

To study this in detail, we again choose the point (yU, M2) = (—500 GeV, 200 GeV). Let 
us consider first the scenario when mj = 380 GeV and m/jo = 110 GeV, so sleptons and Higgs 
bosons are too massive to be on-shell. In this case, the branching ratios are 
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and 

BR(x^ ^ ggX?) = 68%, BRixt ^ q'qxl) = 66%, 

BRixl^m) =11%, BRixt ^ li^xl) =34%, (21) 

BR(x^ ^ uuxl) = 21%. 

The cascade decays dominate, and unfortunately, the on-shell W and Z diagrams result in 
a predominance of hadronic decay products, making it more difficult to identify the primary 
vertex jets. However, when there are only two quark jets, we know that these came from the 
primary vertex, and of the 2333 events that pass the ^/^ and Oacop cuts, there are still 478 
2-jet events. Applying the same analysis to these events as was applied to the cascade events 
of region 2b in Sec. 0, we find that the squark masses can be determined to 4.8 GeV at 
95% CL. We therefore find that despite the predominance of cascades with hadronic decay 
products, we again expect to measure the squark masses to an accuracy of about 1%. 

It is tempting to try to improve our statistics by using some of the events with more 
than 2 jets. For these events, we can try to reconstruct the on-shell Ws and Zs from jet 
pair invariant masses and thereby determine which quark jets are produced at the primary 
vertex. With six jets and the detector resolution of 50% assumed above, it is very difficult to 
determine with any certainty which two of the fifteen possible quark pairs have the correct 
invariant masses to be W and Z decay products. However, for the 996 four jet events, the 
determination is much easier. For these events, we accept only those events where exactly 
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one of the six quark pairs has invariant mass within 10% of either mw or mz- 29% of the 4 
jet events pass this invariant mass cut, and of these, fewer than 1% have misidentified l^s or 
Zs. Therefore, in 29% of the and all of the 2q events we can identify the primary quarks. 
Of these 779 events, 67% are double cascade events. We then calculate m™"^ for each event 
and assume m^^^^^ — Tn^±. Demanding that this hypothesis be kinematically consistent 
leaves 647 events, of which 80% are double cascades. The In C fit to the m™"^ distribution 
gives the squark mass to 4.4 GeV at 95% CL, a shght improvement of our previous result. 

It is important to note that backgrounds to the 4-jet events are large and may make the 
4-jet events difficult to isolate. In particular, if it is necessary to cut events with a pair of 
jets whose invariant mass is near the W or Z mass to reduce the background, the analysis 
of the preceding paragraph is of course not possible. In addition, if /i° is light enough to 
be on-shell in the cascade diagrams, b quarks will dominate the decay products and one 
cannot anti-tag bottom quarks. However, by using a combination of 6-tagging and mass 
reconstruction one might hope to isolate the primary vertex quarks. On the other hand, the 
presence of light sleptons would make the analysis easier, since there are then more lepton 
decay products and more 2-jet events. 

Measurements in region 2b with high and high ruho will be slightly degraded for the 
same reasons as in 2a, namely many cascades with few leptonic decays. One would again 
like to put the large fraction of events with more than two jets to use. Unfortunately, for 
these events, one will generally not be able to identify the primary quark jets, as there 
are no on-shell W and Z decays to reconstruct. However, it may still be possible to get 
information from these events. For example, noting that jets from cascade decays have a 
smooth energy distribution concentrated at lower energies, if one plots the distribution of 1 
jet energies, one might be able to discern the endpoints of the primary decay jet distribution 
above the tail from the cascade jets. In the pessimistic high m^- case we are considering, most 
events will include many jets, and a realistic study would require an accurate simulation of 
hadronization effects and jet reconstruction. 

VIII. ADDITIONAL REMARKS AND CONCLUSIONS 

In our whole analysis above, we have set tan (3 = 2. Raising tan f3 has little effect on the 
squark branching ratios in the /j, > part of the parameter plane, and its effect on branching 
ratios in < regions is to make then more similar to > regions. For example, when 
tan P = 20, the only large effect is that the branching ratio to the LSP in region 3 grows to 
~ 80%. Thus, the decay patterns are nothing new, and by exploring all regions of the plane 
for tan = 2, we have simultaneously roughly analyzed the case of higher tan f3. 

In this paper we have begun to explore the prospects for measuring squark masses at 
future e+e" coUiders, within the context of the MSSM. We have shown that, even if squarks 
choose complex decay patterns, these machines do offer opportunities for making squark 
mass measurements at the level of a few GeV. Such precision measurements would be in- 
valuable for probing the mechanism of supersymmetry breaking in deeper underlying theo- 
ries. Interesting questions for future research include how more realistic physics simulations 
would affect this picture and how squark masses would be studied in the case that gluinos 
are the dominant decay channel. 
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FIGURES 



FIG. 1. The number of squark pairs of the first two generations produced at a 500 GeV e+e~ 
colUder with polarized beams and an integrated luminosity of 10 fb~^ for each beam polarization. 
The four helicity combinations plotted are j^e^ — > qL,RQ- 

FIG. 2. The {iJ,,M2) plane divided into regions with similar squark decay channels, given 
tan [3 = 2 and rriq = 220 GcV. The squark decays of the various regions are described in the 
text. The points of parameter space that we consider in detail in the text are marked. The 
condition ruq < M3 is true only above the dotted line at M2 = 67 GeV. 

FIG. 3. The jet energy distribution resulting from squark decays in region 1 at the point 
(jU, M2) = (—500 GeV, 300 GeV). Each jet produced in such events is binned individually. Cuts 

and detector resolution effects alter the flat shape and make determination of the endpoints more 
difficult. The solid (dashed) histogram represents events with (e^) polarized beams. The 
integrated luminosity assumed is 5 fb^^ per polarization, and the bin size is 4 GeV. 

FIG. 4. Determination of minimum kinematically-allowed squark mass from the two visible 
quark momenta and the LSP mass. The momenta label the particles of q{p\) Xi{P3)'l{P4) ^^'^ 
q{P2) — Xi{P5)q{P6)- The momenta of the two undetected LSPs are constrained to lie on the circle 
C. 

FIG. 5. The distribution of m™"\ the minimum allowed squark mass for a given event, in region 
1 at the point (//, M2) = (—500 GeV, 300 GeV). The distribution for (e^) polarized beams is 
given by the solid (dashed) histogram and is sharply peaked at the actual qi (qR) mass of 220 (210) 
GeV. The integrated luminosity assumed is 5 fb"^ per polarization, and the bin size is 5 GeV. 

FIG. 6. Contours of the branching ratios a) BR(ul uxi) and b) BR(ur uxi) in percent 
in the (/i, M2) plane. In the gaugino-like region 2, ul decays predominantly via cascades, while ur 
is seen to decay primarily to Xi even though other on-shell decays are kinematically allowed. 

FIG. 7. Jet energy distributions at the point (/i, M2) = ( — 500 GeV, 200 GeV) in region 2a. 
Every jet produced in the Monte-Carlo squark decay events is binned individually. The three 
components are a large peak of soft quarks emanating from squarks decaying to X2 ^^'^ xf (long 
dashes) , a flat distribution from direct LSP decays (solid) , and a wide hump from hadronic W and 
Z decays (short dashes). The total distribution is given by the dotted histogram. The analysis is 
independent of the soft jets. The integrated luminosity assumed is 5 fb~^ per polarization, and 
the bin size is 4 GeV. 
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FIG. 8. The total jet energy distribution at the point (/U,M2) = (-500 GeV, 200 GeV) (sohd), 
and the bottom quark energy distribution resulting from the top quark pair production background 
(dashes). In the bottom quark spectrum, we have removed all b quarks produced in top events 
in which at least one b quark has been successfully tagged. We assume mtop = 150 GeV and a 
6-tagging efficiency of 80% |^]. 

FIG. 9. The Feynman diagrams for three-body X2 xf decays mediated by h^, charged 
sleptons and sneutrinos, squarks, and Z and W bosons. 

FIG. 10. Jet energy distributions at the point (/U,M2) = (-100 GeV, 700 GeV) in region 3. 
The solid (dashed) histogram represents events with e~[ {e~^) polarized beams. The soft hump of 
secondary vertex jets is well-separated from the primary vertex jet distribution. The integrated 
luminosity assumed is 5 fb~^ per polarization, and the bin size is 4 GeV. 

FIG. 11. The (;U,M2) plane divided into regions with similar squark decay channels, given 
tan/? = 2 and rriq = 400 GeV. The squark decays of the various regions are as in Fig. I and are 
given in the text. The point in region 2a that we consider in detail is marked. The condition 
niq < M3 is true only above the dotted line at M2 = 121 GeV. 
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